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ABSTRACT 

In  this  paper  we  present  a  novel  digital  filter  intended  for 
quadrature  sub-sample  delay  estimation  based  on  the  Hilbert 
(i complex )  filtering  and  fractional  delaying  of  an  incoming  real- 
valued  signal.  It  is  a  Hilbert  Transform  Filter  (HTF)  having 
variable  fractional  delay  (VFD).  The  filter  is  based  on  a  pair  of 
rotated  VFD  filters  in  the  Farrow  structure.  It  is  capable  of 
performing  the  Hilbertian  as  well  as  VFD  filtering  of  the  incoming 
discrete-time  signal  at  the  same  time.  Thus  one  can  substitute  a 
hitherto  used  cascade  of  the  F1TF  and  the  VFD  filters  with  an 
aggregated  HTF-VFD  filter  proposed  here.  The  technique  is 
simple  to  implement.  The  advantages  lie  in  lower  total  delay 
introduced  by  the  compound  filter  and  in  modular  structure,  which 
is  composed  of  the  same  linear-phase  rotated  VFD  Farrow  sub- 
filters.  The  overall  rotated  VFD  filters  in  the  pair  differ  only  in  the 
value  of  one  parameter  -  the  VFD.  The  proposed  filter  can  be 
applied  to  adaptive  quadrature  sub-sample  estimation  of  delay. 

I.  Introduction 

Fractional  delay  fdters  (FDFs)  find  a  variety  of 
applications  in,  e.g.,  sound  synthesis,  timing  adjustment  in 
digital  receivers,  image  zooming  and  many  others,  which 
have  been  widely  described  so  far  see,  e.g.,  [12].  Recently  a 
FDF  having  variable  fractional  delay  (VFD)  has  been 
applied  to  adaptive  sub-sample  delay  estimation  using  a 
quadrature  detector  based  on  the  FTilbert  (complex)  filtering 
and  fractional  delaying  of  an  incoming  real-valued  signal 
[1],  [2],  [3].  In  the  abovementioned  references  both  of  these 
operations  have  been  realized  by  two  independent  filters, 
namely  the  Hilbert  transform  filter  (HTF)  and  the  VFD 
filter.  In  this  paper  we  propose  to  substitute  these  filters 
with  a  single  aggregated  one,  further  called  the  HTF-VFD, 
having  the  same  functionality,  but  reduced  transport  delay. 
For  this  aim  we  use  polynomial  interpolation.  The 
polynomial  interpolation  can  be  implemented  efficiently 
using  the  idea  based  on  the  Farrow  structure  [4],  [5].  This 
structure  is  very  attractive  because  it  can  provide  variable 
signal  delay  by  changing  only  one  parameter.  This  makes 
on-line  delaying  straightforward  and  feasible.  The  signal 
passes  through  FIR  Farrow  subfilters  either  symmetric  or 
anti-symmetric  see  [4],  hence  their  phase  response  linearity, 
and  is  multiplied  by  appropriate  powers  of  delay  to  produce 
the  output. 

Organization  of  the  paper  is  the  following.  Section  II 
provides  indispensable  background  information  about  the 
fractional  delay  filter  and  the  Farrow  structure.  Section  III 


formulates  the  procedure  of  designing  the  HTF  having 
fractional  delay.  Section  IV  provides  illustrative  design 
examples.  Section  V  summarizes  the  results. 

II.  Background 


Fractional  Delay  Filter 


The  frequency  response  of  an  ideal  FDF  is  defined  as 
Dx  (eJW )  =  exp(-y'K)T ),  |oi|  <  7t  (1) 

where  co  =  2nFT  is  the  normalized  angular  frequency  in 
rad  per  sample,  F  stands  for  the  physical  frequency  in  Hz 
and  T  stands  for  the  sample  interval  in  seconds.  The 
impulse  response  of  this  filter  is 


dx  [«] = 


t~—  \Dx(e>m)eimndm  = 
2n  J 

-n 


simi(«-x) 
n  (n  -t) 


=  sinc(«  -x ) 


(2) 


In  the  case  of  an  FIR  approximation  of  length  N,  x  is 
defined  as  a  sum  of  a  transport  delay  of  the  system  and 
introduced  fractional  delay  a 

x=(N- l)/2  +  a  (3) 


where  the  values  of  a  are  restricted  to  the  interval 
-0.5  <  a  <0.5  in  order  to  deal  with  central  interpolation  - 
the  most  accurate  one. 

The  transfer  function  of  an  FIR  FDF  in  a  direct  form  is 
defined  as 


N-\ 

DNa(z)  =  YjdNaMz" 


n= 0 


(4) 


where  dNa  [«] ,  n  =  0, 1, . . . ,  N  - 1  stands  for  the  impulse 
response  of  the  FIR  filter.  This  filter,  in  case  of  variable 
value  of  a  ,  is  called  the  VFD. 


Farrow  Structure 


High  efficiency  of  a  VFD  can  be  achieved  by  applying 
the  Farrow  structure  [4].  The  main  advantages  of  this 
structure  is  that  the  coefficients  of  its  subfilters  are  either 
symmetrical  or  anti-symmetrical  and  fixed,  independent  of 
the  value  of  introduced  fractional  delay  a  see  [4].  It  means 
that  we  can  change  the  fractional  delay  a  without 
redesigning  the  filter. 
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In  order  to  apply  the  Farrow  structure  we  have  to 
express  the  equation  (4)  in  a  form  of  a  polynomial  of 
fractional  delay  a  .  The  coefficients  of  the  filter  can  be 
rewritten  as 


M 

dNo.[n\  =  'YjCdnl&k  >  n  =  0,l,...,N-l  (5) 

k=0 

where  M+ 1  stands  for  a  number  of  sub  filters  in  the  Farrow 
structure. 

From  the  above  equations  we  have 


M  N- 1  M 

DNa  O)  =  ^^ct[n]z“"a  k  =  ^Q.O)a 1 


where 


k=0  n= 0 


N- 1 


k=0 


Ck(z)  =  ^Jck{n)z  k  =  0,1,...,M 


(6) 


(7) 


n= 0 


are  the  transfer  functions  of  subfilters,  which  in  the  original 
Farrow  work  [4]  are  implemented  in  direct  form.  This 
structure  gives  us  a  flexibility  to  change  the  delay  by 
manipulating  only  one  coefficient,  the  aforementioned 
fractional  delay  (FD). 

x[n 


Fig.  1.  VFD  in  the  Farrow  structure;  y[n\  =  x(n  —  (N  —  l)/2  -a) 
and  d  is  given  by  (10b).  The  subfilters  are  linear-phase  [4], 


III.  Filter  design  and  implementation 


Hilbert  Transform  Filter 


The  frequency  response  of  an  ideal  Hilbert  transform 
filter  having  the  so-called  generalized  linear  phase-response 
is  defined  [6]  as 


H^)  = 


J2exp(-y(© -tt /2)P),  a>e(0,ji) 

0,  toe  (— Tt: ,  0) 


(8) 


where  P  stands  for  the  delay  introduced  by  the  filter.  We 


can  write  the  impulse  response  of  this  filter  using  a  pair  of 
VFDs  [6] 


J  (-l)"/2£/p  /2[«  /2],  77  =  0,  ±  2,  +  4,  .  . 


Therefore  the  design  of  an  FIR  approximation  of  length  2 N 
for  this  ideal  HTF  requires  two  VFDs  of  length  N  each. 

A  general  scheme  of  implementing  this  HTF  having 
variable  fractional  delay  as  an  FIR  filter  in  the  Farrow 
structure  is  presented  in  Fig.  2.  The  total  delay  of  this  filter 
expressed  in  terms  of  a  is 

P=(V-l  +  2a  (10a) 


and  the  fractional  delay  to  be  varied  is 

d  =  2a -1/2  (10b) 

In  Fig.  2  in  the  lower  branch  we  need  to  introduce  an 
additional  delay  z_1  necessary  to  interlace  the  coefficients 
of  the  filter. 

In  case  of  Lagrangian  polynomial  approximation  the 
coefficients  of  a  maximally  flat  VFDcan  be  obtained  using 
the  following  expression 


JV— 1 


‘  Na 


w-n 


k=0 

k*n 


a  -k 
n  -k 


(11) 


for  7?  —0,1,...,  N  1 ,  where  a=([3-/V  +  l)/2  from  (10a). 


x[n] 


a  -1/2 


Fig.  2.  Block  scheme  for  implementing  the  HTF-VFD  in  one  stage 
with  variable  fractional  delay  d  =  2a -1/2  (10b). 

The  procedure  of  designing  the  filter  in  Fig.  2  is  as 
follows.  From  equation  (11)  we  compute  the  coefficients  of 
the  VFD  in  direct  form  and  then,  using  equations  (4),  (6) 

and  (7),  we  compute  the  coefficients  1  .  The  next 

step  is  multiplication  of  these  coefficients  by  appropriate 
powers  of  (-1) ,  as  given  by 

{?A-['C=0  =t-1)"ct[«C  (12) 

in  order  to  get  their  “rotated”  version  (cf.  (9)).  The  resulting 
filter  is  called  in  Fig.  2  the  VFDR  -  variable  fractional 
delayer  “rotated”.  Noteworthy  is  that  we  introduce  the  delay 
a  into  the  upper  branch  in  Fig.  2  and  a  -1/2  into  the 
lower  branch.  The  output  of  the  HTF-VFD  is  a  complex 
signal,  y[n\ ,  described  as 

y[ri\  =  Re  y[ri\  +  j  Imy[??]  (13) 

The  transfer  function  of  the  resulting  HTF-VFD  having 
variable  fractional  delay  2a  is 

JV— 1  JV— 1 

^HTF-VFDW  =  X^(z2)a"  +>_1E^(z2)(a 

k= 0  k= 0 


IV.  Design  examples 

In  this  Section  we  present  some  examples  of  designed 
HTF-VFD  filters. 

A.  Maximally  flat  HTF-VFD 

Example  1.  The  coefficients  of  Farrow  sub  filters  with  a 
maximally  flat  VFDR  of  length  N  =2  for  the  HTF-VFD 
from  Fig.  2  of  length  4,  computed  in  accordance  with  the 
description  as  above  are  the  following 
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{coM/Lo  =  {l/2,  - 1/2}  and  {ct[n]\'n=()  =  {- 1,  - 1}  . 
Consequently 

C0(z)  =  1/2 -l/2z_1  and  Q(z)  =  -l-z^1 , 
and 

^HTF-VFD  (z)  =  ^  Q-  (z  )a  +  Jz  ^  Q-  (z  )(a  _  t) 

A-=0  A=0  Z 

=  1/2 -  a  +  y(l -ex)z_1  - (l/2  +  a)z~2  -y'az~3 
The  impulse  response  and  the  coefficients  of  7/htf-vfd(z) 
for  an  arbitrary  value  of  N  for  this  aggregated  filter  and  also 
for  the  corresponding  cascade  can  be  readily  obtained  using 
the  symbolic  function  HTF-VFD_symb  (N)  written  in 
MATLAB  and  given  below,  where  small  “d”  plays  the  role 
of  “a  ”  in  (14).  An  exemplary  display  of  results  for  N=  2  is 
also  provided. 

function  [htr,htrcascade]=HTF -VFD  symb(N); 

% 

%  program  function  for  the  generation  of  the  coefficients  of 
%  FIR  fractional  delaying  Hilbert  transfonn  filter 
%  HTF-VFD  maximally  flat  around  omega=pi/2 
%  using  symbolic  expressions  for  the  impulse  response 
%  coefficients  htr(n) 

%  as  a  function  of  FD  d;  it  is  assumed  that  -0.5=<d<=0.5 
%  input 

%  N  -  impulse  response  length  of  the  VFD  solely 
%  output 

%  htr  -  the  aggregated  HTF-VFD  coefficients  vector  of  d 

%  htrcascade  -  the  HTF  and  VFD  cascade  imp.  resp.  of  d 

%  use  [htr,htrcascade]=  HTF-VFD_symb(2);  to  start  with  N=2 

syms  d  hi  h2  hi  1  h21  htr  hfsd  hfsdl  hfsd2  htotal 

syms  htrafo  hfsd  sum  hd  g 

%  HTF-VFD  aggregated 

hl=fdfsymb(N);  h2=fdfsymb_shifted(N); 

forn=l:N;  hi l(n)=((-l)A(n-l))*hl(n);  h21(n)=j*((-l)A(n- 

l))*h2(n);  end 

htr(l:2*N)=0;  forn=l:N;  htr(2*n-l)=hll(n);  htr(2*n)=h21(n);  end 

%Cascade  of  linear-phase  HTF  and  VFD 

%Linear-phase  HTF 

hupper=vfd(N,0.5 );  hlower=vfd(N,0); 

for  n=l:N;  hi  lupper(n)=((-l)A(n-l))*hupper(n); 

h2 1  lower(n)=j*((-l)A(n-l))*hlower(n);  end 

h0(l:2*N)=0;  forn=l:N;  h0(2*n-l)=hl lupper(n); 

h0(2*n)=h211ower(n);  end 

htrafo(l  :2*N-l)=hO(2:2*N); 

hfsdl=fdfsymb(N);  %VFD 

for  n=l:2*N-l;  hfsd2(n)=0;  end; 

forn=l:N;  hfsd2(2*n-l)=hfsdl(n);  end 

for  n=l:2:2*N-l;  hfsd(n)=imag(((-l)A(n/2)))*hfsd2(n);  end; 

%Cascade 

%  convolution  of  htrafo  and  hfsd  is  computed  symbolically 
Nx=2*N-l;  hd=htrafo;  g=hfsd;  sum(l)=hd(l)*g(l); 
forn=2:Nx;  sum(n)=hd(l)*g(n); 

for  k=2:n;  sum(n)=hd(k)*g(n-k+l)+sum(n);  end;  end 
sum(2*Nx-l)=hd(Nx)*g(Nx); 

for  n=2*Nx-2:-l:Nx+l;  p=l;  sum(n)=hd(Nx)*g(Nx-p); 

fork=Nx:-l:n-(Nx-l);  sum(n)=hd(k)*g(n-k+l)+sum(n);  end; 
p=p+l ;  end;  htrcascade=sum; 


HTF  -  VFD  imp_resp=conj  (htr') 

HTF  -  VFD_cascade_imp_resp=conj  (htrcascade') 

function  h=fdfsymb(N) 

syms  d  h  hi;  tr=(N-l)/2;  for  n=l:N;  h(n)=l; 

fork=l:N;  if  k~=n  h(n)=h(n)*(tr+d-(k-l))/(n-k);  end;  end;  end 

function  h=fdfsymb_shifted(N) 

syms  d  h  hi;  tr=(N-l)/2;%transport  delay 

for  n=l:N;  h(n)=l; 

fork=l:N;  if  k~=n  h(n)=h(n)*(tr+d-l/2-(k-l))/(n-k);  end;  end;  end 

function  h=vfd(N,de) 
taug=(N-l)/2+de;  forn=l:N;  h(n)=l; 

fork=l:N;  if  k  ~=n h(n)=h(n)*(taug-(k-l))/(n-k);  end;  end;  end 

The  display  of  results  for  N=  2 
HTF-VFD  imp  resp  = 

[  1/2-d];  [  i*(-d+l)]; 

[  -1/2-d];  [  -i*d]; 

HTF-VFD_cascade  imp  resp  = 

[  l/2*i*(l/2-d)];  [  -1/2+d]; 

[  -l/2*i*(l/2-d)+l/2*i*(- 1/2-d)];  [  1/2+d]; 

[  -l/2*i*(-l/2-d)]; 

Magnitude  Response  in  Linear  Scale 


Fig.  3  presents  the  performance  of  the  designed 
narrowband  aggregated  filter,  where  d  is  defined  in  (10b). 
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Fractional  Group  Delay  Response 


Complex  Approximation  Error  Magnitude  in  Pass  Band 


Fig.  3.  The  performance  of  an  aggregated  HTF-VFD  of  length  4, 
maximally  flat  at  f  =  n/  2  ,  from  Example  1. 

For  phase  and  group  delay  responses  in  Fig.  3  the 
transport  delay  having  the  value  of  1.5  sample  intervals  is 
removed  for  better  visualization  of  the  FD.  The  FD  values 
are  set  here  from  the  region  a  e<  0.  0.5  >  for  the  upper 

branch  of  Fig.  2  and  a  -1/2  e<  -  0.5,  0  >  for  the  lower 
branch,  with  an  increment  equal  to  0.1.  It  finally  gives  the 
VFD  value  d  =  2a  -1/2  e<  -0.5,  0.5  >  as  depicted  in  Fig. 
3.  It  is  worth  noticing  that  magnitude  responses  are 
displayed  in  the  range  /e<-l/2, 1/2),  where 
/  =  co  /(2tc  ) ,  on  the  normalized  frequency  axis  and  other 
responses  are  displayed  within  the  range  /  e<  0, 1/2) .  The 
latter  interval  covers  the  pass  band  of  the  filter,  where  these 
characteristics  are  the  object  of  our  interests  and  neglects 
the  stop  band,  where  only  the  attenuation  is  important. 

For  the  sake  of  comparison  Fig.  4  presents  the 
performance  of  the  corresponding  cascade  filter  see  [1],  [2] 
or  [3],  for  the  same  FD  delay  values  as  in  Example  1.  The 
target  filter  length  is  5.  The  transport  delay  introduced  by 
the  cascade  is  2  sample  intervals,  thus  is  greater  than  in 
Example  1.  Also  the  width  of  pass  band  is  unfavorably 
narrower  than  for  the  aggregated  design  proposed  here. 
Thus  in  this  example  we  have  shown  that  the  novel  1TTF- 
VFD  filter  from  Fig.  3  performs  better  than  its  hitherto 
cascade  counterpart  from  Fig.  4,  both  based  on  the  same 
VFD  module.  This  holds  independently  of  the  VFD  module 
length. 

B.  Minimax  F1TF-VFD  optimized  in  two-stage  design 

Finally,  in  Fig.  5,  we  present  the  performance  of  a 
wideband  FITF-VFD  with  a  pair  of  minimax  VFDs 
optimized  in  two-stage  design  [7],  [8].  Two-stage  design 


means  that  at  the  front  end  of  processing  a  linear-phase 
half-band  filter  (FIBF)  produces  two-times  over-sampled 
signal  which  is  passed  through  a  half-band  VFD.  The  result 
is  finally  down-sampled  back  to  the  base-band  frequency. 
In  this  arrangement  the  requirements  for  the  VFD  are 
relaxed  relative  to  one-stage  approach  allowing  for  a 
smaller  number  of  taps  in  the  Farrow  structure.  In  [7]  it  was 
shown  that  the  two-rate  structure  could  be  realized 
efficiently  as  two-stage  one-rate  with  two  parallel  branches. 

Magnitude  Response  in  Linear  Scale 


Normalized  Frequency  1 


Fig.  4.  The  perfonnance  of  a  cascade  of  maximally  flat  F1TF 
and  VFD. 

The  specification  on  VFDs  given  in  [7]  is  the  following. 
The  bandwidth  is  0.45  expressed  in  terms  of  the  normalized 
frequency  / .  The  maximal  absolute  value  of  the  complex 
approximation  error  (CAE)  magnitude  defined  as  the 
difference  between  the  ideal  frequency  response  (1)  and  its 
minimax  approximate  is  0.0042  which  is  equivalent  to 
-47.5  dB.  The  coefficients  of  Farrow  VFD  subfilters 
rounded  to  four  decimal  digits  are  given  in  Figure  4  in  [7]. 
The  number  of  these  subfilters  as  well  as  their  length  is 
jointly  optimized  [7].  The  FIBF  is  a  Nyquist  filter  [11]  of 
order  58,  optimal  in  the  Chebyshev  sense,  with  pass  band 
and  stop  band  edge  normalized  frequencies:  0.45  and  0.55, 
respectively,  requiring  15  taps  of  different  values.  The 
overall  VFD  needs  14  taps  in  parallel  branches.  The 
transport  delay  introduced  by  this  wideband  filter  is  6.5 
sample  intervals. 
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Magnitude  Response  in  Linear  Scale 


Magnitude  Response  Zoomed  in  Pass  Band 


Fig.  5.  The  performance  of  an  aggregated  minimax  HTF-VFD  of 
length  14  with  VFDs  optimized  in  two-stage  design. 

The  maximal  value  of  the  CAE  magnitude  of  the 
resulting  HTF-VFD  (defined  as  the  absolute  value  of  the 
difference  between  the  ideal  frequency  response  (8)  and  its 
approximate)  arranged  as  shown  in  Fig.  2  is  approximately 
-51.5  dB  in  the  bandwidth  of  0.225  located  symmetrically 


around  the  normalized  frequency  0.25.  Note  that  in  case  of 
d—  0,  the  target  HTF-VFD  filter  performance  can  be 
assessed  in  a  way  shown  in  [9],  [10]  or  [11]  on  the  basis  of 
the  performance  of  the  VFD  used  in  Fig.  2  as  the  VFDR. 

V.  Conclusions 

A  novel  efficient  design  of  a  digital  filter  for  quadrature 
sub-sample  delay  estimation  was  proposed  and  described  in 
this  paper.  The  coefficients  of  this  filter  are  obtained  from 
the  coefficients  of  a  pair  of  identical  VFDs  appropriately 
“rotated”  and  interlaced.  Maximally  flat  (Lagrangian)  [6] 
approximation  as  well  as  minimax  approximation  [7],  [8]  of 
VFDs  optimized  in  two-stage  design  was  employed  to 
compute  coefficients  of  the  designed  HTF-VFD  filter. 
Using  in  this  design  the  structure  invented  by  Farrow  [4] 
gives  a  possibility  of  changing  the  value  of  introduced  delay 
in  real-time,  without  redesigning  the  coefficients.  The 
original  Farrow  structure  [4]  is  composed  of  linear-phase 
subfilters  and  as  such  can  be  used  straightforwardly  without 
any  need  for  modifications.  (This  seems  unnoticed  by  many 
authors.)  It  was  shown  that  the  proposed  aggregated  HTF- 
VFD  exhibits  wider  bandwidth  and  introduces  smaller 
integer  delay  in  comparison  with  a  hitherto  applied  cascade 
design  based  on  a  series  connection  of  two  independent 
filters:  the  linear-phase  HTF  and  the  VFD. 

References 

[1] D.L.  Maskell  and  G.S.  Woods,  “Adaptive  subsample  delay  estimation 
using  a  modified  quadrature  phase  detector”,  IEEE  Trans.  Circuits  Syst.  - 
II,  vol.  52,  No.  10,  Oct.  2005,  pp.  669-674. 

[2] D.L.  Maskell,  G.S.  Woods,  and  A.  Kerans,  “A  hardware  efficient 
implementation  of  an  adaptive  subsample  delay  estimator”,  in  Proc. 
ISCAS  ’2004,  vol.  Ill,  pp.  3 1 7-320 

[3] D.L.  Maskell  and  G.S.  Woods,  “Adaptive  subsample  delay  estimator 
using  quadrature  estimator”,  Electr.  Letters,  vol.  40,  No.  5,  4th  March  2004, 
pp.  347-349. 

[4] C.W.  Farrow,  “A  continuously  variable  delay  element”,  in  Proc.  IEEE 
Int.  Symp.  Circuits  Syst.,  vol.  3,  Espoo,  Finland,  June  7-9,  1988,  pp.  2641- 
2645. 

[5] S.  Samadi,  A.  M.  Omair  and  M.N.S.  Swamy,  “Results  on  maximally  flat 
fractional-delay  systems”,  IEEE  Trans.  Circuits  Syst.  -  I,  vol.  51,  No.  11, 
Nov.  2004,  pp.  2271-2286. 

[6] E.  Hermanowicz,  Special  Discrete-Time  Filters  and  Applications. 
EXIT,  Warsaw,  Poland,  2005,  Chapter  4. 

[7] E.  Hermanowicz  and  H.  Johansson,  “On  designing  minimax  adjustable 
wideband  fractional  delay  FIR  filters  using  two-rate  approach”,  in  Proc. 
European  Conference  on  Circuit  Theory  and  Design,  ECCTD  2005,  Cork, 
Ireland,  29  August  -  1  September  2005,  vol.  I,  pp.  441-444. 

[8] H.  Johansson  and  E.  Hermanowicz,  “Adjustable  fractional-delay  filters 
utilizing  the  Farrow  structure  and  multirate  techniques”,  in  Proc.  of  The 
Sixth  International  Workshop  on  Spectral  Methods  and  Multirate  Signal 
Processing,  SMMSP2006,  Florence,  Italy,  2-3  September  2006. 

[9] P.P.  Vaidyanathan  and  T.Q.  Nguyen,  “A  trick  for  the  design  of  FIR 
half-band  filters”,  IEEE  Trans. Circuits  Syst.  vol.  CAS-34,  Mar.  1987,  pp. 
297-300. 

[10]  P.P.  Vaidyanathan,  “Design  and  implementation  of  digital  FIR  filters”, 
in  Handbook  on  Digital  Signal  Processing,  Edited  by  D.F.  Elliott, 
Academic  Press  Inc.,  1987,  pp.  140-143. 

[1 1] P.P.  Vaidyanathan,  Multirate  Systems  and  Filter  Banks,  Prentice  Hall, 
1993,  pp. 185-186. 

[12]  T.  I.  Laakso,  V.  Valimaki,  M.  Karjalainen  and  U.  K.  Laine:  Splitting 
the  unit  delay.  Tools  for  fractional  delay  filter  design.  IEEE  Signal 
Processing  Magazine,  January  1996,  pp.  29-60. 


©2007  EURASIP 


1057 


EUSIPCO,  Poznan  2007 


